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SUMMARY 


This paper presents a procedure to allow temperature dependent mechanical properties 
to be used in the Engineering Analysis Language (EAL) System for solid structural 
elements. The implementation of this procedure is accomplished by including a modular 
runstream in the main EAL runstream. The procedure is applicable for models with 
multiple materials and with anisotropic properties, and can easily be incorporated into an 
existing EAL runstream. The procedure is applicable for EAL elastic solid elements and is 
described in detail in the paper, followed by a description of the validation of the 
procedure. A listing of the EAL runstream used to validate the procedure is included in the 
Appendix. 


INTRODUCTION 

Finite element analysis is widely used to study various thermal and mechanical 
problems. One finite element program that is frequently used is the Engineering Analysis 
Language (EAL) System (ref. 1). EAL is a modular system of individual analysis 
processors which may be used in any appropriate sequence to perform a variety of 
analyses. In performing thermal analyses, EAL allows the use of temperature dependent 
thermal properties, i.e., thermal conductivity, specific heat, density, and emissivity. 
However, in the solution of thermal stress problems, EAL limits the use of temperature 
dependent mechanical properties to one- and two-dimensional isotropic materials. For 
orthotropic or anisotropic three-dimensional materials, EAL does not allow for the 
inclusion of temperature dependent mechanical properties, e.g., modulus of elasticity, 
shear modulus, coefficient of thermal expansion, and Poisson’s ratio. In many thermal 
stress problems, mechanical properties may vary greatly in space and/or time, and thus it is 
necessary to include this temperature dependence to insure accurate results. 
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The purpose of this paper is to present a procedure to implement temperature 
dependent mechanical properties for solid structural elements in EAL. The implementation 
of this procedure is accomplished by including a few modifications to an existing EAL 
runstream. The procedure is applicable for models with multiple materials and with 
anisotropic properties, and is written for EAL S81, S61, and S41 elements, but can be 
extended to other elements. After a brief description of EAL, the implementation procedure 
is described, followed by a description of the procedure. A listing of the EAL runstream 
used to validate the procedure is included in the Appendix. 

SYMBOLS 

DNS8 
E 

E(x) 

E43 
EAL 
Ell 

FLUX K81 
L 

PROP BTAB 
RECI 
S41 
S61 
S81 
T(x) 

T 0 


name given to data set of temperature dependent densities 
EAL processor used in structural analysis 
Young’s modulus 
plate elements 

Engineering Analysis Language 
EAL processor 

EAL table generated by processor TAFP containing element 

temperatures 

rod length 

EAL name for table containing flexibility matrix 

EAL funct on to calculate reciprocal 

4-node tetrahedron element 

6-node wedge element 

8-node isoparametric brick element 

temperature 

rod temperature at x = 0 

EAL processor that calculates average element temperatures 


TAFP 
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TK81 name given to data set of average element temperatures 

u(x) displacement 

x coordinate in rod length direction 

XNT1 interpolation function in EAL 

* EAL symbol preceding all executive control statements 

$ signifies comment statement in EAL 

! signifies register action command (used for simple mathematics) in EAL 

Greek symbols 
a stress 

P coefficient in linear variation of elastic modulus with temperature 

y coefficient in linear variation of thermal expansion with temperature 

a linear coefficient of thermal expansion 

DESCRIPTION OF EAL 

The Engineering Analysis Language (EAL) System is a high-level finite-element 
program which evolved from the SPAR computer program (ref. 2). The EAL System 
contains individual processor; which communicate through a data base consisting of one or 
more libraries of data sets. A sequence of processor execution commands, which are 
written using an executive control language, is referred to as a runstream. Data sets 
typically contain data describing the finite-element model (e.g., material properties and 
nodal coordinates), as well as response information such as displacements, stresses, and 
temperatures. The EAL system uses flexible, FORTRAN-like statements which allow 
branching, looping, testing of data, and calling of runstreams (similar to calling a 
FORTRAN subroutine). The procedure presented here is a modular runstream which is 
called from within a main EAL runstream (similar to a FORTRAN main program with 
subroutines). 
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PROCEDURE FOR IMPLEMENTATION OF 
TEMPERATURE-DEPENDENT PROPERTIES 

The runstream used to determine temperature dependent mechanical properties for 
solid structural elements and to calculate the structural response for two example problems 
is listed in the Appendix. The procedure for implementing the temperature dependent 
properties includes a runstrt am which creates tables of each property as a function of 
temperature for each material (called DATA TABL in the Appendix) and an additional 
runstream to calculate the temperature dependent properties for each element (called TMDP 
MAPR in the Appendix). In general, the runstream TMDP MAPR defines element 
temperatures (average of the element nodal temperatures), interpolates individual flexibility 
terms for element temperatures, and replaces individual element stiffness matrices stored in 
the data base. Implementing the procedure in EAL requires modifications to an existing 
EAL runstream to provide tables of properties as a function of temperature and a call to a 
newly created runstream (TMDP MAPR) which calculates the temperature dependent 
properties. Details of the procedure required to calculate the temperature dependent 
properties are discussed in this section. 

To use temperature dependent mechanical properties, one must first provide a 
property table called PROP BTAB for each different material used (see runstream DATA 
TABL in the Appendix). These are the tables that are normally used to input constant 
mechanical properties, and no change is required in the structure of these tables. However, 
it is important that the density of each material be unique, since it is the density that is used 
in the present runstreams to differentiate between materials with different properties in this 
procedure. Hence, if two different materials happen to have the same density, it will be 
necessary to use “dummy densities” in the PROP BTAB tables. Since the density values 
from the PROP BTAB tables are not used in the structural analysis, there is no risk of 
using incorrect data. The actual density is obtained from a separate table of temperature 
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dependent densities. After generating the PROP BTAB tables for each material, a unique 
table of temperature dependent properties is created for each of the non-zero terms. Thus, 
up to 31 temperature dependent tables may be required for each material since there are 31 
entries in PROP BTAB. Although the modulus may be linear, EAL stores the flexibility 
coefficients which are related to the inverse of the moduli. Hence, it is imperative that 
enough data points be included so that a linear interpolation between two points produces a 
satisfactory approximation to the inverse of the modulus. 

After generating the tables of temperature dependent properties (tables labeled CDEN, 
CA11, CA21,..., etc.), the runstream TMDP MAPR is called to insert the temperature 
dependent properties into the EAL data base between the execution of the processors E and 
EKS in the structural analysis. A listing of the runstream TMDP MAPR is included in the 
Appendix, and the steps followed are outlined in the flow chart in figure 1. In the 
description that follows, the term vector is used to denote a one-dimensional list of values. 

The first step in the runstream TMDP MAPR is to obtain the element temperatures. 
This can be done in one of two methods. If a thermal analysis has been performed 
previously using EAL (and thus a STAT TEMP data set exists), the processor TAFP can be 
executed to calculate the element temperatures. Once the element temperatures are obtained 
with processor TAFP, a separate vector of element temperatures (called TK81 in the 
Appendix) must be created. This vector must be created if processor TAFP is used since 
the element temperatures that are created by processor TAFP are part of another table (i.e., 
FLUX K81 if K81 brick elements are used). (The method using TAFP is included in the 
runstream in the Appendix, however, it is commented out and not used since a thermal 
analysis was not performed prior to the structural analysis.) If a thermal analysis does not 
precede the structural analysis, the element temperatures can be input directly into a vector 
called TK81, or nodal temperatures can be input and averaged to obtain the element 
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temperature vector TK81. The procedure used in the runstream listed in the Appendix 
assumes no prior thermal analysis has been performed. 

The remaining portion of runstream TMDP MAPR should be transparent to the user, 
and is now discussed. If elements other than S81, S61 or S41 are used, or more than two 
materials are used, this portion of runstream TMDP MAPR would need to be modified. 
With the calculation of the element temperatures complete, the element temperatures are 
now used to interpolate the temperature dependent properties for each element. The EAL 
function XNT1 is used in TMDP MAPR to linearly interpolate between the element 
temperature (stored in TK81) and the tables with temperature dependent properties. This is 
done for every property for each material and results in a vector for each temperature 
dependent property for each material, assuming all the same type of elements are of the 
same material. Thus, if two materials are used with S81 elements, there will be a vector for 
each property (up to 31) assuming all S81 elements are material 1, and a second set of 
vectors assuming all S81 elements are material 2. 

It is necessary to operate on the temperature dependent property vectors previously 
created to obtain one vector for each property for each element type. This vector will 
contain values for a given property as a function of temperature and material. To create this 
vector, the material dependent densities of each element type are first extracted from the 
EAL data base. This is done using the Ell processor. After the densities are extracted, 
they need to be placed in a single vector (called DNS8 in the Appendix) containing the 
constant, but material dependent, densities. To illustrate this, consider a density vector 
extracted with the Ell processor with 3 elements and 2 different materials represented by 
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DNS8 ■ 


Pi 

P2 

Pi 


The density vector DNS8 is now used to obtain unit vectors for each material. The unit 
vector for material 2 will have a 1 for each element with material 2 and a 0 for each element 
of any other material. The unit vector for material 1 will have a 1 for each element with 
material 1 and a 0 for each element of any other material. To obtain the unit vector for 
material 2, subtract pj from each density in DNS 8 as follows 


“ “ ™ 




mm ^ 

Pi 


Pi 


0 

P2 

- 

Pi 

= 

P2 ' Pi 

Pi 


Pi 


0 


Then, using the EAL function RECI, multiply each entry by its reciprocal, noting that in 
EAL, the reciprocal of 0 is 0. 


mm 


r i 

0 


t 

0 

P2 ‘ Pi 


1/(P2-Pl) 

0 


0 



__ _ 


0 

1 

0 


This is the unit vector for material 2. It has a 1 for each element of material 2, and a 0 for 
each element of material 1 . To obtain the unit vector for material 1 , subtract the unit vector 
for material 2 from the identity vector, as shown. 
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1 


0 


1 

1 

- 

1 

= 

0 

1 


0 


1 


A similar procedure would follow if more that two materials are used. 


These unit vectors are now used to obtain a vector of temperature dependent 
properties for each element This is done by multiplying the unit vector for each material 
by the property vector that was generated by assuming all of that element type were of that 
material. To illustrate this, consider the coefficient of thermal expansion in the x direction 
as a function of temperature assuming all S81 elements are made of material 1. The 
coefficient of thermal expansion vector calculated previously for material 1 is of the form 

AL1C = o,,= [ <x xl (T) a xl (T) a„i(T) ] 

Likewise, assuming all the elements are of material 2, the coefficient of thermal expansion 
vector is given by, 

AL1W s a x2 = [ a x2 (T) a x2 (T) a x2 (T) ] 

Now multiply the material 1 coefficient of thermal expansion vector by the material 1 unit 
vector to redefine AL1C as 


AL1C ^ [ a xl (T) a xl (T) a xl (T) ] 


0 

1 


[ 0 a x i (T) 0 ] 


0 


and multiply the material 2 coefficient of thermal expansion vector by the material 2 unit 
vector to redefine AL1 W as 
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AL1W = [ a x2 (T) a x2 (T) a x2 (T) ] 


0 


= [ a x2 (T) 0 a x2 (T) ] 


Now sum the individual coefficient of thermal expansion vectors as follows 
AL18 * [ 0 a xl (T) 0] + [ ot x2 (T) 0 a x2 (T) ] 

= [ a x2 (T) a x ,(T) a x2 (T) ] 

to obtain a vector of the temperature dependent coefficient of thermal expansion for all S81 
elements. This same procedure is repeated for each of the non-zero properties in the 
flexibility matrix. 


These property vectors are now put into a table (called NEW PROP in the Appendix) 
with 31 columns (the number of entries in the PROP BTAB table) and the number of rows 
equal to the number of elements of a given type. Entries that are zero need not be put in 
this table. This temperature dependent flexibility matrix is now embedded into the EAL 
data set by the Ell processor. Now all of the S81 elements have the correct temperature 
dependent mechanical properties. If S61 or S41 elements are used in the runstream, the 
same procedure should be repeated for each element type. If elements other than S81, S61, 
or S41 are used, and temperature dependent properties are required, a similar procedure 
could be followed. However, the manner in which the properties are saved (and thus 
extracted and embedded) within EAL in the E state is different. 


EXAMPLE PROBLEMS USED FOR PROCEDURE VALIDATION 

To validate the procedure using temperature dependent mechanical properties, two 
example problems were considered. One example compares stresses in a thin plate 
modeled using E43 (two-dimensional plate) elements with those which assume a three- 
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dimensional model consisting of S81 (three-dimensional solid) elements. A second 
example compares stresses in a rod calculated using S81 elements with a closed-form 
analytical solution. In both examples, the S81 elements used temperature dependent 
mechanical properties calculated using the present procedure. 

Rat Plate Constrained Along All Four Edges and Subject to a Bi-linear 

Temperature Distribution 

This validation compares stresses in a thin plate modeled with S81 solid elements 
with the same plate modeled with E43 plate elements. The plate is unconstrained in the 
through the thickness direction and is constrained in both the x and y in-plane directions. 
An illustration of the plate, including element numbers and nodal temperatures, is shown in 
figure 2. Two different materials are used in the analysis for generality. The reason for 
comparing an E43 element solution with an S81 element solution is that EAL allows 
temperature dependent isotropic properties for E43 elements. Thus, since the plate was 
modeled assuming isotropic S81 elements, one can compare the current temperature 
dependent mechanical properties solution using S81 elements with an established EAL 
solution using temperature dependent E43 elements. Although this does not allow the 
validation of all of the entries in the flexibility matrix (since EAL only allows isotropic 
temperature dependent E43 elements), it does allow the validation of the present procedure 
for isotropic properties. Since the procedure is the same for anisotropic properties, it is 
assumed that if the procedure works correctly for isotropic properties, the procedure also 
works correcdy for anisotropic properties. 

The temperature dependent properties used for the thin plate are given in Table 1, and 
the resulting stresses from the S81 and E43 solutions are compared in Table 2. The 
comparison of the stress results shown in Table 2 indicate that stresses calculated using 
S81 solid elements (temperature dependent properties obtained using the present procedure) 
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compare quite favorably with those calculated using E43 plate elements (temperature 
dependence of the properties calculated internally in EAL). A comparison of stresses, 
where the temperature dependant properties (for S81 elements) are interpolated between 
entries 1000 °F apart as well as 100 °F apart, is also shown in Table 2. Young’s modulus 
for the E43 elements is obtained directly in the EAL runstream, rather than from the inverse 
of the modulus as with S81 elements. Thus the E43 element solution should not be as 
sensitive to property interpolation as long as the properties vary nearly linearly with 
temperature. It is anticipated that the 100 °F interpolation would be more accurate for the 
S81 solution than the 1000 °F interpolation. However, since the stresses are compared 
only with another EAL solution, nothing can be said about the accuracy of the solution, 
only that the solutions from the two methods are consistent. 

Rod Constrained at Both Ends and Subject to a Linear Temperature Distribution 

The second method of procedure validation involves comparing stresses from a rod 
modeled using S81 elements with stresses from a corresponding closed-form analytical 
solution. The temperature distribution in the rod, which is constrained at both ends, is 
assumed to be a linear function of x as 

T(x) = T 0 + J x I (1) 

where x is the coordinate along the length of the rod, Tq is the temperature at x = 0, AT is 
the temperature difference between the ends of the rod, and L is the length of the rod. In 
addition, the modulus of elasticity E and coefficient of thermal expansion a are assumed to 
be linear functions of temperature. Since the temperature is a linear function of x, both E 
and a are linear functions of x, given as 

E(x) = Eq [1 + p T(x)] = Eo[l + P |t 0 + *-Pj] 


( 2 ) 



<x(x) = aoYT(x) = oto y [t 0 + 
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(3) 


where P and y are constants. Using the symbolic mathematics program Mathematica (ref. 
3), an analytical closed form solution is obtained for the displacement in the rod and is 
given as 


Lcc 0 Y 


u(x) = 


[ T » 


x AT' 


3 AT 


L a 0 (3 T 0 2 + 3 T 0 AT + AT 2 ) y log 1 + P T 0 + 

+ 3 [ log[ 1 + p T 0 ] - log[l + P (T 0 + AT)]] 

+ { L a 0 Y [ T 0 3 log (l + P T 0 ) + 3 T 0 3 AT log (l + p T 0 ) 


+ 3 T 0 AT 2 log (l + P T 0 )+AT 3 log(l + P T 0 ) 


- T 0 3 log (l +p(T 0 + AT)) ]) 

/ {3 AT [log (1 + P T 0 ) - log (1 + P (T 0 + AT)) ]) (4) 

Though the above equation for the displacement is not valid for AT = 0, it does approach 
the correct displacement in the limit as AT — > 0. From the displacement, the stress in the 
rod is calculated from 

o(x) = E(x) - E W "W T « (5) 

to be 


E 0 a 0 p y AT (3 T 0 2 + 3 T 0 AT + AT 2 ) 

3 [log (1 + P T 0 ) - log (1 + P (T 0 + AT)) ] 


( 6 ) 


Equation (6) shows that the stress is constant for the linear temperature, modulus, and 
coefficient of thermal expansion assumed in this case. Equation (6) is not valid at AT = 0, 


13 


i.e., for a constant temperature rod. However, in the limit as AT — » 0, eq. (6) approaches 
the correct stress. The constant terms used in eqs. (4, 6) are 


Eq — 1.0 x K^psi 

(7a) 

1.0 x 10' 6 in. 

(7b) 

a 0 " in. °F 

T 0 = 200.0 °F 

(7c) 

AT = 800.0 °F 

(7d) 

1000 °F 

(7e) 

1 

y_ 1000°F 

(70 

L = 1 in. 

(7g) 


The analytical displacement obtained from eq. (4) is compared with the EAL solution using 
20 elements in Table 3. The mechanical properties in the EAL solution are interpolated 
from property data every 100 °F. Although the percent error between the analytical and 
EAL displacements is not zero everywhere, it is small, with the largest error being 4.15 
percent. The analytical stress obtained from eq. (6) is compared with the EAL stress using 
20 elements in Table 4. As can be seen in Table 4, when the rod is assumed to be a 
constant temperature of 200 J F or 1000 °F, eq. (6) and EAL give exactly the same stress. 
When the temperature distribution in the rod is linear, the EAL solution depends on the 
accuracy of interpolating the inverse of the modulus. If only the end point values are used 
for interpolation (i.e., 200 °F and 1000 °F), the EAL stress is -619.8 psi, compared with 
the exact value of -647.3 psi. However, when the table for the inverse of the modulus 
includes values at 100 °F intervals, the EAL stress is -646.7 psi. This increased accuracy is 
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due to the fact that the inverse of the modulus is nearly linear between values in the table 
spaced at 100 °F intervals. 


CONCLUDING REMARKS 

A procedure to allow the Engineering Analysis Language (EAL) System to use 
temperature dependent mechznical properties has been presented and shown to be accurate. 
The procedure involves calling a run stream (subroutine) from the main EAL runstream to 
calculate the temperature dependent properties. Though the procedure described here has 
been written for the EAL S81, S61, and S41 elements, it can be modified for other 
elements. A modular approach has been taken in the development of the procedure, thus 
enabling an existing EAL runstream to be easily modified. 
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APPENDIX 


This Appendix contains a listing of the EAL runstream used to calculate the stresses 
in the validation of the procedure assuming temperature dependent mechanical properties. 
The stresses can be calculated in either a rod or a plate with this runstream. The stresses in 
the plate can be calculated using either E43 plate elements or S81 isoparametric “brick” 
elements. The procedure used to calculate the temperature dependent mechanical properties 
is contained in a subroutine called TMDP MAPR. In the listing, two different materials are 
used, and are differentiated by “C” and “W”. Comments are preceded by the symbol, $, 
and are generally in lower case. 


DD'NN 

TITLE' TEST MODEL 

This is a test program. It will calculate the temperature dependent 
' thermal stresses in a plate using either S81 elements or E43 
’ elements. Change the parameter GEOM to get either plate (E43) or 
' solid (S81) elements. Since the E43 elements only allow temperature 
' dependent isotropic properties, temperature dependent isotropic 
' properties are also used for the S81 elements. 

$ 

♦ECHO 1,2, 3, 4 
•ABORT 1 
•XQTU1 
*CM= 150000 


^******************************* 4 ^*************************************** 


$ This subroutine does the structural analysis 

$ 

J * * * ** * * * 4* ** * * * 4- * ** * * * * * * * * * + * * * * * * * * * * 4» * * * * * * * * * * + * * * * * * * + * + + * Hi * * * * * * * * 


$ 


• (29 STRU ANAL) ZSTRU 

$ 


•XQTTAB 

$ 


1GEOM = 0 $ 0 for plates, 1 for solids 

IDIMN = 2 $ 1 or 2 depending on dimension of problem desired 

$ Determine temperature for TEMP NODE data set 


$ 


1TINC = 40.0 

$ temperature increment for each node (2-D problem) 

ITND1 = 40.0 

$ temperature of node 1 (2-D problem) 

!TX0 = 200.0 

$ temperature at one side of rod (1-D problem) 

!TXL = 1000.0 

$ temperature at other side of rod (1-D problem) 

!XLEN= 1.0 

$ length in x direction 

!YLEN= 1.0 

$ length in y direction 

IZLEN = 0.1 

$ thickness 

!NX = 6 

$ number of nodes in x direction 

!NY = 7 

$ number of nodes in y direction (odd number) 

!NZ = 2 

$ number of nodes in z direction 
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!NELX = NX - 1 
!NELY = NY - 1 
!NELZ = NZ - 1 
!NNOD = NZ * NX * NY 
START "NNOD" 4 56 


$ number of elements in x direction 
$ number of elements in y direction 
$ number of elements in z direction 
$ total number of of nodes 


INI = 1 

!N2 = NX * NY + 1 
!DX » XLEN/NELX 
!DY = YLEN/NELY 
!DZ = ZLEN/NELZ 
1XINC = 1 

$ Generate the nodes 
JLOC 

FORMAT=l 

"Nl" 0.0.0. "XLEN” 0. 0. "NX" "XINC" "NY" 

"NX" 0. "YLEN" 0. "XLEN" "YLEN" 0. 

"N2" 0. 0. "DZ" "XLEN" 0. "DZ" "NX" "XINC" "NY" 
"NX” 0. "YLEN” "DZ" "XLEN" "YLEN" "DZ" 


$ 

* DC ALL (29 STRU BCON) 

$ 

MATC $ for plate elements 

1 1.0E+6 0.3 0.69 1.0E-6 

2 1.0E+6 0.3 0.69 1.0E-6 
SA 

1 "ZLEN" 

NMAT=2 

2 "ZLEN" 

$ 

♦XQTAUS 

* DC ALL (29 DATA TABL) $ insert properties as a function of temperature 

S 

♦XQTDCU 
PRINT 1 MATC 
PRINT 1 NODA TEMP 
$ 

♦XQTELD 
1IINC = 1 
1JINC = NX 
1KINC = NX * NY 
!NELY = NY- 1/2 
!N2 = NX * NELY + 1 
*IF("GEOM" EQ 0): *JUMP600 
S81 

NSECT = 1 


GROUP = 1 

1 "NELX" "NELY" "NELZ" "IINC" "JINC" "KINC" 0 1 
NSECT = 2 

"N2 M "NELX" "NELY" "NELZ" "IINC" "JINC" "KINC" 0 1 
♦JUMP 610 
♦LABEL 600 
$ 

E43 

NSECT=1 


!J1 = 1 

!J2 = J1 +IINC 
!J3 = J2 + JINC 
!J4 = J1 + JINC 



"J1" "J2" ”J3" "J4" 1 "NELX" "NELY" 

NSECT=2 

!J1 = NX * NELY + 1 
!J2 = J1 + IINC 
!J3 = J2 + JINC 
!J4 = J1 + JINC 

"Jl" "J2” "J3" "J4" 1 "NELX" "NELY" 

$ 

♦LABEL 610 
$ 

♦XQTU1 

*TI (1 TD MATC 1 1) 

0.0 1.494E+7 0.4 2.00E-6 

1000.0 1.300E+7 0.3 3.26E-6 

2000.0 1.146E+7 0.2 5.28E-6 

3000.0 1.025E+7 0.1 6.31E-6 

*TI (1 TD MATC 2 1) 

0.0 5.903E+7 1.054E-1 1.00E-6 

1000.0 3.712E+7 1.034E-1 2.26E-6 

2000.0 2.685E+7 1.046E-1 3.28E-6 

3000.0 2.103E+7 8.433E-2 4.31E-6 

*TI (1 TD THICKNESS 1 1) 

"ZLEN" 

"ZLEN" 

$ 

♦XQT DCU 

*IF("GEOM" EQ 1): PRINT 1 DEF S81 
$ 

♦XQTE 

T=l.E-20 -.IE-2 0.01 0.01 40. 0.01 0.01 0.01 
$ 

♦IFfGEOM" EQ 1): * DC ALL (29 TMDP MAPR) $ calculate temp. dep. material prop. 
$ 

♦XQTEKS 

•IFC’GEOM” EQ 1): *JUMP 720 
RESET TLIB=1 
RESET N3T=1 
♦LABEL 720 
SRESET SRFZ=0.0 
$ 

♦XQT SEQ 
$ 

♦XQT TAN 
$ 

♦XQTK 
*ONLINE=2 
♦XQT DRSI 
♦ONLINE=l 
$ 

♦XQTEQNF 
$ 

♦XQT SSOL 
$ 

♦XQT VPRT 
PRINT STAT DISP 1 1 
PRINT STAT REAC 1 1 
$ 


$ to avoid imminent singularity error 
$ joint elimination sequence 


$ equivalent nodal forces 
$ displacements, reactions 
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$ 

♦XQTPSF 

DIV 1000. 1000. 1000. 1000. 

$ 

•XQTES 

$ 

♦XQTAUS 

DEFI S=STAT DISP 1 

GDSP=LTOG(S) 

$ 

•XQT DCU 

CHANGE 1 GDSP AUS 1 1 PATR DISP 1 1 
PRINT 1 JLOC BTAB 
TOC 1 
$ 

♦RETURN 

* ZSTRU 


$ 

$ 


j************************************************************************ 


$ 

$ This subroutine calculates temperature dependent material properties 
$ 


j************************************************************************ 


$ 

* (29 TMDP MAPR) ZTMDP 

$ 

$ Create a vector of element center temperatures and then interpolate 

$ to find the corresponding material properties 

•XQTAUS 

!NK81 = TOC NJ(1 S81 EFIL MASK MASK) 

!NK81 = 


INELM = NX - 1 
INELM = NY - 1 * NELM 
!NELM = 


!NY = 

!NX = 

TABLE(NI= 1 NJ= "NK8 1 ") : TK81 1 

*IF("DIMN" EQ 2): *JUMP5132 $ vary temperatures in 1 direction 
IDELT = NY - 1.0 
IDELT = TXL - TX0 / DELT 
INLY = NY - 1 
IINK = 1 
•LABEL 7 196 

*IF("INK" GT "NLY"): *JUMP 7197 

111 = INK - 1 * NELX + 1 

112 = INK * NELX 

1TELM = INK - 1.0 * DELT + TX0 
ITELN = INK • DELT + TX0 
ITELM = TELM + TELN / 2.0 
J="I1","I2": "TELM" 

IINK = INK + 1 


•JUMP 7196 
♦LABEL 7197 
•JUMP 5133 
$ 

$ 

•LABEL 5132 
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’"LABEL 900 

♦IFfND" GT "NK81"): *JUMP 910 
!I13 = DS,13,"ND", 1(1, DEF,S81, 18,8) 

!I14 = DS,14,"ND”, 1(1, DEF.S81, 18,8) 

!I15 = DS,15,"NDM(1,DEF,S81,18,8) 

!I16 = DS,16,"ND",1(1,DEF,S81,18,8) 

!TPND = 113 + 114 + 115 + 116 /4 * TINC $ average element temperature 

*IF("TPND" GT 0.1): J="ND": "TPND" 

♦IF("TPND" LT 0. 1 ): J=''ND" : "TND 1" $ for constant temperature (i.e., TINC=0) 

!ND = ND + 1 
•JUMP 900 
♦LABEL 910 
♦LABEL 5133 
$ 

$ 

$ If one calculates nodal temperatures in EAL, use the following procedure to get element temperatures 
$*XQT TAFP $ to get element center temperatures 

$ TEMP=STAT TEMP 
$ COMPUTE ELEMENT FLUXES: K81 
$ 

$*XQT DCU 
$ PRINT 1 FLUX K81 
$ 

$*XQT AUS 

$ 1NK81 = TOC NJ(1 S81 EFIL MASK MASK) 

$ DEFI F1=FLUX K81 1 1 
$ TABLE(NI=1 J NJ=“NK81”): TK81 
$ TRAN(SOUR=Fl ,ILIM=ULIM=“NK81”,SSKI=2) 

$ 

$ Now the element temperatures (TK81) have been embedded from one of the 2 methods, print them out 
♦XQTDCU 
PRINT 1 TK81 
$ 

$Now use those element temperatures to interpolate and find properties of S81's 
♦XQT AUS 

DENC=XNT1(CDEN,TK81) $ interpolate to find the density 

DENW=XNT1 (WDEN.TK8 1) 

A11C=XNT1(CA11,TK81) $ interpolate to find all 

A1 1W=XNT1(WA1 1.TK81) 

A21C=XNT1(CA21,TK81) $ interpolate to find a21 

A21W=XNT1(WA21,TK81) 

A22C=XNT 1 (C A22.TK8 1 ) $ interpolate to find a22 

A22W=XNT1(WA22,TK81) 

A31C=XNT1(CA31 ,TK81) $ interpolate to find a31 

A3 1 W=XNT1(WA3 1 ,TK81 ) 

A32C=XNT 1 (C A32.TK8 1 ) $ interpolate to find a32 

A32W=XNT 1 (WA32.TK8 1 ) 

A33C=XNT 1 (C A33 ,TK8 1 ) $ interpolate to find a33 

A33W=XNT1(WA33,TK81) 

A44C=XNT1(CA44,TK81) $ interpolate to find a44 

A44W=XNT1(WA44,TK81) 

A55C=XNT1(CA55,TK81) $ interpolate to find a55 

A55W=XNT1(WA55,TK81) 

A66C=XNT1(CA66,TK81) $ interpolate to find a66 

A66W=XNT1 (WA66.TK8 1 ) 

AL1C=XNT1(CAL1,TK81) $ interpolate to find alphal 
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AL1 W=XNT1(WAL1 ,TK8 1 ) 

AL2C=XNT1 (C AL2.TK8 1 ) $ interpolate to find alpha2 

AL2W=XNT1(WAL2,TK8 1) 

AL3C=XNT1(CAL3,TK81) $ interpolate to find alpha3 
AL3 W=XNT 1 (WAL3 ,TK8 1 ) 

YXXC=XNT1(CYXX,TK81) $ interpolate to find Yxx 
YXXW=XNT1(WYXX,TK81) 

YYYC=XNT1 (CYYY.TK8 1) $ interpolate to find Yyy 

YYYW=XNT1 (WYYY.TK8 1) 

YZZC=XNT1(CYZZ,TK81) $ interpolate to find Yzz 
YZZW=XNT1(WYZZ,TK81) 

YXYC=XNT 1 (C YXY.TK8 1 ) $ interpolate to find Yxy 

YXYW=XNT 1 (WYX Y,TK8 1 ) 

YXZC=XNT1(CYXZ,TK81) $ interpolate to find Yxz 
YXZW=XNT 1 (WYXZ.TK8 1 ) 

YYZC=XNT1(CYYZ,TK81) $ interpolate to find Yzy 
YYZW=XNT1(WYYZ,TK81) 

$ 

*XQT Ell $ Get the vector of material densities 

EXTRACT 

SOURCE = S81 

CONT SPEC: SECTION 1 $ density only 

CREATE DNS8 S81 1 1 $ vector containing densities of S81's 

$ 

♦XQTDCU 

PRINT DNS8 S81 1 1 $ print the density vector 

$ 

$ Set up unit vectors to be used in the calculation 
*XQT AUS 

DEFI DN=DNS8 S81 1 1 $ material dependent densities 

TABLE(NI= 1 ,NJ="NK8 1 "): CARB VECT: J=1,"NK81": "CRHO" 

DEFI CV=CARB VECT 
VI = SUM(DN,-1.0 CV) 

DI = RECI(V1) 

VT = PROD(V 1 ,DI) $ unit vector 

TABLE(NI=U^J="NK81"): UNIT VECT: J=l,"NK81 n : 1.0 
DEFI UV=UNIT VECT 
VC = SUM(VT,-1.0 UV) 

VC = PROD(VC,- 1 .0 UV) $ unit vector 

$ v ™=" - ■ - ------- 

$ Now multiply the unit vectors by the temperature dependent properties and sum to get temperature and 
$ material dependent properties 
*XQT AUS 

DENC = PROD(DENC.VC) 

DENW = PROD(DENW,VT) 

DEN8 = SUM(DENC,DENW) $ vector of temperature dependent densities 

A1 1C = PROD(Al 1C, VC) 

A1 1W = PROD(Al 1W.VT) 

A1 18 = SUM(A1 1C.A1 1 W) $ vector of temperature dependent al 1 

A21C = PROD(A21C,VC) 

A21W = PROD(A2 1 W,VT) 

A218 = SUM(A21C,A21W) $ vector of temperauire dependent a21 

A22C = PROD(A22C,VC) 

A22W = PROD(A22W,VT) 

A228 * SUM(A22C,A22W) $ vector of temperature dependent a22 

A31C = PROD(A31C,VC) 

A31W = PROD( A3 1 W, VT) 

A318 = SUM(A3 1C, A3 1W) $ vector of temperauire dependent a31 



A32C = PROD(A32C,VC) 
A32W = PROD(A32W,VT) 
A328 = SUM(A32C,A32W) 
A33C = PROD(A33C,VC) 
A33W = PROD(A33W,VT) 
A338 = SUM(A33C,A33W) 
A44C = PROD(A44C,VC) 
A44W = PROD(A44W,VT) 
A448 = SUM(A44C,A44W) 
A55C = PROD(A55C,VC) 
A55W = PROD(A55W,VT) 
A558 = SUM(A55C,A55W) 
A66C = PROD(A66C,VC) 
A66W = PROD(A66W,VT) 
A668 = SUM(A66C,A66W) 
AL1C = PROD(ALIC.VC) 
AL1W = PROD(AL 1 W,VT) 
AL18 = SUM(AL1C,AL1W) 
AL2C = PROD(AL2C,VC) 
AL2W = PROD(AL2W,VT) 
AL28 = SUM(AL2C AL2W) 
AL3C = PROD(AL3C,VC) 
AL3W = PR0D(AL3W,VT) 
AL38 = SUM(AL3C,AL3W) 
YXXC = PROD(YXXC,VC) 
YXXW = PROD(YXXW,VT) 
YXX8 = SUM(YXXC,YXXW) 
YYYC = PROD(YYYC,VC) 
YYYW = PROD(YYYW,VT) 
YYY8 = SUM(YY Y C, YY YW) 
YZZC = PROD(YZZC,VC) 
YZZW = PROD(YZZW,VT) 
YZZ8 = SUM(YZZC.YZZW) 
YXYC = PROD(YXY C,V C) 
YXYW = PROD(YXYW,VT) 
YXY8 = SUM(YXYC,YXYW) 
YXZC = PROD(YXZC,VC) 
YXZW = PROD(YXZW,VT) 
YXZ8 = SUM(YXZC.YXZW) 
YYZC = PROD(YYZC,VC) 
YYZW = PROD(YYZW,VT) 
YYZ8 = SUM(YYZC,YYZW) 
$ 


$ vector of temperature dependent a32 
$ vector of temperature dependent a33 
$ vector of temperature dependent a44 
$ vector of temperature dependent a55 
$ vector of temperature dependent a66 
$ vector of temperature dependent alphal 
$ vector of temperature dependent alpha2 
$ vector of temperature dependent alpha3 
$ Yxx 
$ Yyy 
$ Yzz 
$ Yxy 
$ Yxz 
$ Yyz 


$Now insert each of these vectors into a dataset 


TABLE(NI=3 1 ,NJ="NK8 1 ”): NEW PROP 

TRAN(SOUR=DEN8,ILIM=l JLIM="NK8r,SBAS=0,DBAS=0,DSKI=30) 
TRAN(SOUR=Al 18.ILIM=ULIM="NK8r’,SBAS=0,DBAS=1 ,DSKI=30) 
TRAN(SOUR=A2 1 8,ILIM=1 r FLIM="NK8 1 ”,SBAS=0,DB AS=2,DSKI=30) 
TRAN(SOUR= A228,ELIM= 1 JLIM="NK8 1 " ,SB AS=0,DB AS=3,DSKI=30) 
TRAN(SOUR= A3 1 8,ILIM= 1 r fLIM=”NK8 1 " ,SB AS=0,DB AS=4,DSKI=30) 
TRAN(SOUR=A328,ILIM= 1 JLIM="NK8 1 " ,SBAS=0,DBAS=5,DSKI=30) 
TRAN(SOUR= A338,ILIM= 1 JLIM="NK8 1 " ,SB AS=0,DB AS=6 J)SKI=30) 
TRAN(SOUR= A448JLIM= 1 JLIM="NK8 1 ” ,SB AS=0,DB AS=1 0,DSKI=30) 
TRAN(SOUR= A558,ILIM= 1 JLIM="NK8 1 " ,SB AS=0 J)B AS=1 54^SKI=30) 
TRAN(SOUR=A668JLIM=l JLIM="NK81 ",SBAS=0,DBAS=21 ,DSKI=30) 
TRAN(SOUR=AL 18,ILIM= 1 r FLIM="NK8 1”,SB AS=0,DBAS=22,DSKI=30) 
TRAN(SOUR=AL28,ILIM=l JL1M="NK81",SBAS=0,DBAS=23,DSKI=30) 
TRAN(SOUR=AL38,ILIM=l JLIM="NK8r,SBAS=0,DBAS=24J)SKI=30) 



TRAN(S0UR=YXX8,ILIM= 1 JLIM="NK8 1 \SB AS=0,DB AS=25,DSKI=30) 
TRAN(S0UR=YYY8,ILIM= 1 JLIM="NK8 1 ",SB AS=0,DB AS=26,DSKI=30) 
TRAN(S0UR=YZZ8,ILIM=1 JLIM="NK8r , ) SBAS=0,DBAS=27,DSKI=30) 
TR AN(S0UR=YXY8 JLIM= 1 JLIM="NK8 1 " ,SB AS=0,DB AS=28,DSKI=30) 
TRAN(S0UR=YXZ8,ILIM=1 JL M=”NK81 ”,SB AS=0,DB AS=29,DSKI=30) 
TRAN(S0UR=YYZ8,ILIM=1 JLIM="NK81 ”,SBAS=0,DBAS=30,DSKI=30) 

$ 

$Now embed the just created dataset of properties into Eli 

•XQTEI1 

EMBED 

SOURCE = NEW PROP 
CONT SPEC: SECTION 1,31 
DEST = S81 
$ 

♦RETURN 

* ZTMDP 


$ 

$ 

g********************** ******************* ***+*************************** 

$ 

$ This subroutine sets the boundary conditions for the structural analysis 
$ 

JH^********************t ***********★*********+***++*********+*++****+**** 

$ 

♦ (29 STRU BCON) ZBCON 

CON=l 

$ 

111 = 1 


ZERO 2 
•LABEL 200 

•IFCIl” GT "NX"): * JUMP 100 
!I2 = NX * NELY + II 


!I3 = NX • NY + II 
!I4 = NX * NY + 12 
"II" 

” 12 " 

”13” 

"14" 


111 = 11 + 1 
•JUMP 200 
•LABEL 100 


$ 

ZERO 1 2 3 4 5 6: 1 
•RETURN 

* ZBCON 


$ 

$ Define some of the constraints needed 
$ 


g*********************************************************************** 


$ 

$ Create the tables of data 
$ 


g*** I************************************************************ ******** 

$ 


* (29 DATA TABL) ZDATA 

$ Now input the structural data sets 
$ 



$ Set up the NODA TEMP data set 
TABLE(NI=l,NJ="NNOD"): NODA TEMPI 
*IF("DIMN" EQ 2): *JUMP5128 
!DELT = NY - 1.0 
!DELT = TXL - TXO / DELT 
!INK = 1 
♦LABEL 7192 

♦IFflNK" GT "NY"): *JUMP7193 
III = INK - 1 * NX + 1 
!I2 = INK * NX 

!TELM = INK - 1.0 * DELT + TXO 
J="I1","I2": "TELM" $ 1-D problem 
!INK = INK + 1 
*JUMP 7192 
♦LABEL 7193 

*■»*... 

♦LABEL 7194 

•IFflNKl” GT "NY"): ♦JUMP7195 
!I1 = INK - 1 * NX + 1 
!I2 = INK * NX 

1TELM = INK1 - 1.0 * DELT + TXO 
J=”H","I2": "TELM" 

IINK1 = INK1 + 1 
!INK = INK + 1 
♦JUMP 7194 
♦LABEL 7195 
♦JUMP 5129 
$ 

$ 

♦LABEL 5128 
!ND = 1 

!NN2 = NNOD/2 
♦LABEL 210 

♦IF("ND" GT "NN2"): ♦JUMP 200 
1TMPD = ND * TINC 
♦IFC’TMPD” EQ 0.0): 1TMPD = TND1 
!ND = 

!TMPD = 

J="ND": "TMPD" $ 2-D problem 
!ND = ND + 1 
♦JUMP 210 
♦LABEL 200 

*IF("ND" GT "NNOD"): *JUMP 220 
1TMPD = ND - NN2 * TINC 
♦IFC’TMPD" EQ 0.0): 1TMPD = TND1 
J="ND": "TMPD" 

!ND = ND + 1 
♦JUMP 200 
♦LABEL 220 
♦LABEL 5129 
$ 

1CRHO = 0.2 
1WRHO = 0.6 

$ These are tje property tables normally required by EAL. These numbers will be replaced later. 
TABLE(NI=3 1 ,NJ=2): PROP BTAB 2 21 
J=1 $ Material 1 

”CRHO”> 
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2.0E-6> 

0.0000E-9 2.0000E-6> 
O.OOOOE-9 0.0000E-9 2.000E-6> 


0.0 

0.0 

0.0 

1.0E-6> 

0.0 

0.0 

0.0 

0.0 

1.0E-6> 

0.0 

0.0 

0.0 

0.0 

0.0 1 ,0E-6> 


1.0E-6 1.0E-6 1.0E-6> 

1.0 1.0 1.0 1.0 1.0 1.0$ 

s 

J=2 $ Material 2 
"WRHO"> 

2.121E-8> 

-6.364E-9 2.121E-8> 

-6.364E-9 -6.364E-9 2.121E-8> 


0.0 

0.0 

0.0 

5.605E-8> 

0.0 

0.0 

0.0 

0.0 5.605E-8> 

0.0 

0.0 

0.0 

0.0 0.0 5.605E-8> 


2.969E-6 2.969E-6 2.969E-6> 

1.0 1.0 1.0 1.0 1.0 1.0$ 

$ 

$Set up tables for temperature dependent properties of material 1 

T ABLE(NI=2,NJ=2): CDEN: 1=1,2 $ density 

J=l: 0.0 "CRHO" 

J=2: 3000.0 "CRHO" 

TABLE(NI=2,NJ=20): CA1 1: 1=1,2 $ al 1 = 1/Ej 

J=l: 0.0 6.694E-8 

J=2: 100.0 6.782E-8 

J=3: 200.0 6.872E-8 

J=4: 300.0 6.965E-8 

J=5: 400.0 7.060E-8 

J=6: 500.0 7.158E-8 

J=7: 600.0 7.259E-8 

J=8: 700.0 7.363E-8 

J=9: 800.0 7.469E-8 

J=10: m0 7.579E-8 
J=ll: 1000.0 7.694E-8 

J=12: 1 100.0 7.785E-8 

J=13: 1200.0 7.879E-8 

J=14: 1300.0 7.976E-8 

J=15: 1400.0 8.075E-8 

J=16: 1500.0 8.177E-8 

J=17: 1600.0 8.281E-8 

J=18: 1700.0 8.388E-8 

J=19: 2000.00 8.724E-8 

J=20: 3000.00 9.754E-8 

TABLE(NI=2,NJ=20): CA21: h 1,2 $ a21 = -v 12 /E, 

J=l: 0.0 -2.678E-8 

J=2: 100.0 -2.636E-8 
J=3: 200.0 -2.595E-8 
J=4: 300.0 -2.555E-8 

J=5: 400.0 -2.517E-8 

J=6: 500.0 -2.479E-8 

J=7: 600.0 -2.443E-8 

J=8: 700.0 -2.408E-8 
J=9: 800.0 -2.373E-8 

J=10: m.O -2.340E-8 
J=ll: 1000.0 -2.308E-8 

J= 12: 1100.0 -2.236E-8 
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J=13: 1200.0 

-2.168E-8 

J=14: 1300.0 

-2.104E-8 

J=15: 1400.0 

-2.044E-8 

J=16: 1500.0 

-1.987E-8 

J=17: 1600.0 

-1.934E-8 

J=18: 1700.0 

-1.883E-8 

J=19: 2000.00 

-1.745E-8 

J=20: 3000.00 

-9.754E-9 

TABLE(NI=2,NJ=20): CA22: 

<— 1 

II 

O 

O 

6.694E-8 

J=2: 100.0 

6.782E-8 

J=3: 200.0 

6.872E-8 

J=4: 300.0 

6.965E-8 

J=5: 400.0 

7.060E-8 

J=6: 500.0 

7.158E-8 

J=7: 600.0 

7.259E-8 

J=8: 700.0 

7.363E-8 

J=9: 800.0 

7.469E-8 

J=10: 900.0 

7.579E-8 

J=ll: 1000.0 

7.694E-8 

J=12: 1100.0 

7.785E-8 

J=13: 1200.0 

7.879E-8 

J=14: 1300.0 

7.976E-8 

J=15: 1400.0 

8.075E-8 

J=16: 1500.0 

8.177E-8 

J=17: 1600.0 

8.28 IE-8 

J=18: 1700.0 

8.388E-8 

J=19: 2000.00 

8.724E-8 

J=20: 3000.00 

9.754E-8 


T ABLE(N1=2,NJ=20): CA3 1 : 1= 1 ,2 
J=l: 0.0 -2.678E-8 

J=2: 100.0 -2.636E-8 

J=3: 200.0 -2.595E-8 

J=4: 300.0 -2.555E-8 

J=5: 400.0 -2.517E-8 

J=6: 500.0 -2.479E-8 

J=7: 600.0 -2.443E-8 

J=8: 700.0 -2.408E-8 

J=9: 800.0 -2.373E-8 

J=10: 900.0 -2.340E-8 

J=ll: 1000.0 -2.308E-8 

J=12: 1100.0 -2.236E-8 

J=13: 1200.0 -2.168E-8 

J=14: 1300.0 -2.104E-8 

J=15: 1400.0 -2.044E-8 

J=16: 1500.0 -1.987E-8 

J=17: 1600.0 -1.934E-8 

J=18: 1700.0 -1.883E-8 

J=19: 2000.00 -1.745E-8 
J=20: 3000.00 -9.754E-9 
TABLE(NI=2,NJ=20): CA32: 1=1,2 
J=l: 0.0 -2.678E-8 

J=2: 100.0 -2.636E-8 

J=3: 200.0 -2.595E-8 

J=4: 300.0 -2.555E-8 

J=5: 400.0 -2.517E-8 

J=6: 500.0 -2.479E-8 

J=7: 600.0 -2.443E-8 


$ a22 = 1/E 2 


$a31 =-v 13 /E, 


$ a32 = -V23/E2 
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J=8: 700.0 

-2.408E-8 

J=9: 800.0 

-2.373E-8 

J=10: 900.0 

-2.340E-8 

J— 11: 1000.0 

-2.308E-8 

J=12: 1100.0 

-2.236E-8 

J=13: 1200.0 

-2.168E-8 

J=14: 1300.0 

-2.104E-8 

J=15: 1400.0 

-2.044E-8 

J=16: 1500.0 

-1.987E-8 

J=17: 1600.0 

-1.934E-8 

J=18: 1700.0 

-1.883E-8 

J=19: 2000.00 

-1.745E-8 

J=20: 3000.00 

-9.754E-9 

TABLE(NI=2J4J=20): CA33: 

<— 1 

II 

O 

O 

6.694E-8 

J=2: 100.0 

6.782E-8 

J=3: 200.0 

6.872E-8 

J=4: 300.0 

6.965E-8 

J=5: 400.0 

7.060E-8 

J=6; 500.0 

7.158E-8 

J=7: 600.0 

7.259E-8 

J=8: 700.0 

7.363E-8 

J=9: 800.0 

7.469E-8 

J= 1 0: 900.0 

7.579E-8 

J=ll: 1000.0 

7.694E-8 

J=12: 1100.0 

7.785E-8 

J=13: 1200.0 

7.879E-8 

J=14: 1300.0 

7.976E-8 

J=15: 1400.0 

8.075E-8 

J=16: 1500.0 

8.177E-8 

J=17: 1600.0 

8.281E-8 

J=18: 1700.0 

8.388E-8 

J=19: 2000.00 

8.724E-8 

J=20: 3000.00 

9.754E-8 

TABLE(N1=2,NJ=20): CA44: 

T 

O 

b 

1.874E-7 

J=2: 100.0 

1.886E-7 

J=3: 200.0 

1.898E-7 

J=4: 300.0 

1.910E-7 

J=5: 400.0 

1.922E-7 

J=6: 500.0 

1.935E-7 

J=7: 600.0 

1.948E-7 

J=8: 700.0 

1.960E-7 

J=9: 800.0 

1.973E-7 

J=10: 900.0 

1.987E-7 

J=ll: 1000.0 

2.000E-7 

J=12: 1100.0 

2.009E-7 

J=13: 1200.0 

2.018E-7 

J=14: 1300.0 

2.027E-7 

J=15: 1400.0 

2.036E-7 

J=16: 1500.0 

2.046E-7 

J=17: 1600.0 

2.055E-7 

J=18: 1700.0 

2.065E-7 

J=19: 2000.00 

2.094E-7 

J=20: 3000.00 

2.146E-7 


TABLE(NI=2,NJ=20); CA55: 1=1,2 
J=l: 0.0 1.874E-7 

J=2: 100.0 1.886E-7 


$ a33 = 1/E 3 


$ a44 = I/G23 


$ a55 = 1/G 31 
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J=3: 200.0 1.898E-7 

J=4: 300.0 1.910E-7 

J=5: 400.0 1.922E-7 

J=6: 500.0 1.935E-7 

J=7: 600.0 1.948E-7 

J=8: 700.0 1.960E-7 

J=9: 800.0 1.973E-7 

J=10: 900.0 1.987E-7 

J=ll: 1000.0 2.000E-7 

J=12: 1100.0 2.009E-7 

J=13: 1200.0 2.018E-7 

J=14: 1300.0 2.027E-7 

J=15: 1400.0 2.036E-7 

J=16: 1500.0 2.046E-7 

J=17: 1600.0 2.055E-7 

J=18: 1700.0 2.065E-7 

J=19: 2000.00 2.094E-7 

J=20: 3000.00 2.146E-7 

T ABLE(NI=2,NJ =20): CA66: 1=1,2 $a66=l/G 12 

J=l: 0.0 1.874E-7 

J=2: 100.0 1.886E-7 

J=3: 200.0 1.898E-7 

J=4: 300.0 1.910E-7 

J=5: 400.0 1.922E-7 

J=6: 500.0 1.935E-7 

J=7: 600.0 1.948E-7 

J=8: 700.0 1.960E-7 

J=9: 800.0 1.973E-7 

J=10: 900.0 1.987E-7 

J=ll: 1000.0 2.000E-7 

J=12: 1100.0 2.009E-7 

J=13: 1200.0 2.018E-7 

J=14: 1300.0 2.027E-7 

J=15: 1400.0 2.036E-7 

J=16: 1500.0 2.046E-7 

J=17: 1600.0 2.055E-7 

J=18: 1700.0 2.065E-7 

J=19: 2000.00 2.094E-7 
J=20: 3000.00 2.146E-7 

TABLE(NI=2,NJ=4): CALI: 1=1,2 $ alpha 1 


J=l: 

0.0 

2.00E-6 

J=2: 

1000.0 

3.26E-6 

J=3: 

2000.0 

5.28E-6 

J=4: 

3000.0 

6.31E-6 


TABLE(N1=2,NJ=4): CAL2: 1=1,2 $ alpha2 


J=l: 

000.0 

2.00E-6 

J=2: 

1000.0 

3.26E-6 

J=3: 

2000.0 

5.28E-6 

J=4: 

3000.0 

6.31E-6 


TABLE(NI=2,NJ=4): CAL3: I=i,2 $alpha3 


J=l: 

000.0 

2.00E-6 

J=2: 

1000.0 

3.26E-6 

J=3: 

2000.0 

5.28E-6 

J=4: 

3000.0 

6.31E-6 

TABLE(NI=2,NJ=2): CYXX: 

J=l: 

0.0 

1.0 

J=2: 

30000.0 

1.0 


TABLE(NI=2,NJ=2): CYYY: 1=1,2 $ Yyy 
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J=l: 0.0 1.0 

J=2: 30000.0 1.0 

TABLE(NI=2,NJ=2): CYZZ: 1=1,2 $ Yzz 

J=l: 0.0 1.0 

J=2: 30000.0 1.0 

TABLE(NI=2,NJ=2): CYXY: 1=1,2 $ Yxy 

J=l: 0.0 1.0 

J=2: 30000.0 1.0 

TABLE(NI=2,NJ=2): CYXZ: 1=1,2 $ Yxz 

J=l: 0.0 1.0 

J=2: 30000.0 1.0 

TABLE(NI=2,NJ=2): CYYZ: 1=1,2 $Yyz 

J=l: 0.0 1.0 

J=2: 30000.0 1.0 

$ 

$Set up tables few temperature dependent properties of material 2 

TABLE(NI=2,NJ=2): WDEN: 1=1,2 $ density 

J=l: 0.0 "WRHO" 

J=2: 3000.0 "WRHO" 

TABLE(NI=2,NJ=22): WA1 1: 1=1,2 $ al 1 = 1/Ej 


J= 1: 

.0 

1.694E-08 

J= 2: 

100.0 

1.759E-08 

J= 3: 

200.0 

1.830E-08 

J= 4: 

300.0 

1.906E-08 

J= 5: 

400.0 

1.989E-08 

J= 6: 

500.0 

2.080E-08 

J= 7: 

600.0 

2.179E-08 

J= 8: 

700.0 

2.289E-08 

J= 9: 

800.0 

2.410E-08 

J=10: 

900.0 

2.544E-08 

J— 11: 

1000.0 

2.694E-08 

J=12: 

1100.0 

2.771E-08 

J=13: 

1200.0 

2.852E-08 

J=14: 

1300.0 

2.938E-08 

J=15: 

1400.0 

3.029E-08 

J=16: 1500.0 

3.126E-08 

J=17: 1600.0 

3.230E-08 

J=18: 1700.0 

3.341E-08 

J=19: 

1800.0 

3.459E-08 

J=20: 

1900.0 

3.587E-08 

J=21: 

2000.0 

3.724E-08 

J=22: 3000.00 

4.754E-8 


TABLE(NI=2J S JJ=22): WA21: 1=1,2 $a21=-v 12 /E 1 

J= 1: .0 -1.786E-09 

J= 2: 100.0 -1.852E-09 

J= 3: 200.0 -1.924E-09 

J=4: 300.0 -2.001 E-09 

J= 5: 400.0 -2.085E-09 
J= 6: 500.0 -2.176E-09 

J= 7: 600.0 -2.276E-09 
J= 8: 700.0 -2 385E-09 

J= 9: 800.0 -2.505E-09 

J=10: 900.0 -2.638E-09 

J=1I: 1000.0 -2.786E-09 
J=12: 1100.0 -2.867E-09 
J=13; 1200.0 -2.954E-09 
J=14: 1300.0 -3.046E-09 
J=15: 1400.0 -3.144E-09 
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J=16: 1500.0 -3.248E-09 
J=17: 1600.0 -3.360E-09 
J=18: 1700.0 -3.479E-09 
J=19: 1800.0 -3.607E-09 
J=20: 1900.0 -3.745E-09 
J=21: 2000.0 -3.895E-09 
J=22: 3000.00 -4.009E-9 
TABLE(NI=2,NJ=22): WA22: 1=1,2 

J= 1 : .0 1.694E-08 

J= 2: 100.0 1.759E-08 

J= 3: 200.0 1.830E-08 

J= 4: 300.0 1.906E-08 

J= 5: 400.0 1.989E-08 

J= 6 : 500.0 2.080E-08 

J= 7: 600.0 2.179E-08 

J= 8 : 700.0 2.289E-08 

J= 9: 800.0 2.410E-08 

J= 10: 900.0 2.544E-08 

J=ll: 1000.0 2.694E-08 
J=12: 1100.0 2.771E-08 
J=13: 1200.0 2.852E-08 
J=14: 1300.0 2.938E-08 
J=15: 1400.0 3.029E-08 
J=16: 1500.0 3.126E-08 
J=17: 1600.0 3.230E-08 
J=18: 1700.0 3.341E-08 
J=19: 1800.0 3.459E-08 
J=20: 1900.0 3.587E-08 
J=21: 2000.0 3.724E-08 
J=22: 3000.00 4.754E-8 
TABLE(NI=2,NJ=22): WA31: 1-1,2 
J= 1: .0 -1.786E-09 

J= 2: 100.0 -1.852E-09 

J= 3: 200.0 -1.924E-09 

J= 4: 300.0 -2.001E-09 

J= 5: 400.0 -2.085E-09 

J= 6 : 500.0 -2.176E-09 

J= 7: 600.0 -2.276E-09 

J= 8 : 700.0 -2.385E-09 

J= 9: 800.0 -2.505E-09 

J=10: 900.0 -2.638E-09 

J=ll: 1000.0 -2.786E-09 
J=12: 1100.0 -2.867E-09 
J=13: 1200.0 -2.954E-09 
J=14: 1300.0 -3.046E-09 
J=15: 1400.0 -3.144E-09 
J=16: 1500.0 -3.248E-09 
J=17: 1600.0 -3.360E-09 
J=18: 1700.0 -3.479E-09 
J=19: 1800.0 -3.607E-09 
J=20: 1900.0 -3.745E-09 
J=21: 2000.0 -3.895E-09 
J=22: 3000.00 -4.009E-9 


TABLE(NI=2,NJ=22): WA32: 1=1,2 

J= 1: 

.0 - 

1.786E-09 

J= 2: 

100.0 

-1.852E-09 

J= 3: 

200.0 

-1.924E-09 

J= 4: 

300.0 

-2.001E-09 


$ a 22 = I/E 2 


$ a31 = -Vj 3 /Ej 


$ a32 = -V 23 /E 2 
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J= 5: 

400.0 

-2.085E-09 

J= 6: 

500.0 

-2.176E-09 

J= 7: 

600.0 

-2.276E-09 

J= 8: 

700.0 

-2.385E-09 

J= 9: 

800.0 

-2.505E-09 

J=10: 

900.0 

-2.638E-09 

J-ll: 1000.0 

-2.786E-09 

J=12: 1100.0 

-2.867E-09 

J=13: 

1200.0 

-2.954E-09 

J=14: 

1300.0 

-3.046E-09 

J=15: 

1400.0 

-3.144E-09 

J=16: 

1500.0 

-3.248E-09 

J=17: 

1600.0 

-3.360E-09 

J=18: 

1700.0 

-3.479E-09 

J=19: 

1800.0 

-3.607E-09 

J=20: 1900.0 

-3.745E-09 

J=21: 2000.0 

-3.895E-09 

J=22: 3000.00 

-4.009E-9 


TABLE(NI=2,NJ=22): WA33: 1=1,2 $ a33 = 1/E 3 

J= 1: .0 1.694E-08 

J= 2: 100.0 1.759E-08 

J= 3: 200.0 1.830E-08 

J= 4: 300.0 1.906E-08 

J= 5: 400.0 1.989E-08 

J= 6: 500.0 2.080E-08 

J= 7: 600.0 2. 179E-08 

J= 8: 700.0 2.289E-08 

J= 9: 800.0 2.410E-08 

J=10: 900.0 2.544E-08 

J=ll: 1000.0 2.694E-08 
J=12: 1100.0 2.771E-08 
J=13: 1200.0 2.852E-08 
J=14: 1300.0 2.938E-08 
J=15: 1400.0 3.029E-08 
J=16: 1500.0 3.126E-08 
J=17: 1600.0 3.230E-08 
J=18: 1700.0 3.341E-08 
J=19: 1800.0 3.459E-08 
J=20: 1900.0 3.587E-08 
J=21: 2000.0 3.724E-08 
J=22: 3000.00 4.754E-8 

TABLE(NI=2,NJ=22): WA44: 1=1,2 $ a44 = 1/G 23 


J= 1: 

.0 

3.7450E-08 

J= 2: 

100.0 

3.8889E-08 

J= 3: 

200.0 

4.0443E-08 

J= 4: 

300.0 

4.2127E-08 

J= 5: 

400.0 

4.3957E-08 

J= 6: 

500.0 

4.5953E-08 

J= 7: 

600.0 

4.8138E-08 

J= 8: 

700.0 

5.0543E-08 

J= 9: 

800.0 

5.3200E-08 

J=10: 

900.0 

5.6151E-08 

J— 11: 

1000.0 

5.9450E-08 

J=12: 

1100.0 

6.U47E-08 

J=13: 

1200.0 

6.2943E-08 

J=14: 

1300.0 

6.4848E-08 

J=15: 

1400.0 

6.6872E-08 

J=16: 

1500.0 

6.9026E-08 



J=17: 1600.0 7.1324E-08 
J=18: 1700.0 7.3780E-08 
J=19: 1800.0 7.641 IE-08 
J=20: 1900.0 7.9237E-08 
J=21: 2000.0 8.2280E-08 
J=22: 3000.00 1.031E-7 

T ABLE(NI=2,NJ=22): WA55: 1=1,2 $a55=l/G 31 

J= 1: .0 3.7450E-08 

J= 2: 100.0 3.8889E-08 

J= 3: 200.0 4.0443E-08 

J= 4: 300.0 4.2127E-08 

J= 5: 400.0 4.3957E-08 

J= 6: 500.0 4.5953E-08 

J= 7: 600.0 4.8138E-08 

J= 8: 700.0 5.0543E-08 

J= 9: 800.0 5.3200E-08 

J=10: 900.0 5.6151E-08 

Jail: 1000.0 5.9450E-08 
J=12: 1100.0 6.1147E-08 
J=13: 1200.0 6.2943E-08 
J=14: 1300.0 6.4848E-08 
J=15: 1400.0 6.6872E-08 
J=16: 1500.0 6.9026E-08 
J=17: 1600.0 7.1324E-08 
J=18: 1700.0 7.3780E-08 
J=19: 1800.0 7.641 IE-08 
J=20: 1900.0 7.9237E-08 
J=21: 2000.0 8.2280E-08 
J=22: 3000.00 1.031E-7 

TABLE(NI=2,NJ=22): WA66: 1=1,2 $a66=l/G 12 

J= 1: .0 3.7450E-08 

J= 2: 100.0 3.8889E-08 

J= 3: 200.0 4.0443E-08 

J= 4: 300.0 4.2127E-08 

J= 5: 400.0 4.3957E-08 

J= 6: 500.0 4.5953E-08 

J= 7: 600.0 4.8138E-08 

J= 8: 700.0 5.0543E-08 

J= 9: 800.0 5.3200E-08 

J=10: 900.0 5.6151E-08 

J=ll: 1000.0 5.9450E-08 
J=12: 1100.0 6.1147E-08 
J=13: 1200.0 6.2943E-08 
J=14: 1300.0 6.4848E-08 
J=15: 1400.0 6.6872E-08 
J=16: 1500.0 6.9026E-08 
J=17: 1600.0 7.1324E-08 
J=18: 1700.0 7.3780E-08 
J=19: 1800.0 7.641 IE-08 
J=20: 1900.0 7.9237E-08 
J=21: 2000.0 8.2280E-08 
J=22: 3000.00 1.031E-7 

TABLE(NI=2,NJ=4): WAL 1; 1= 1 ,2 $ alpha 1 


Jal: 

0.0 

1.00E-6 

J=2: 

1000.0 

2.26E-6 

J=3: 

2000.0 

3.28E-6 

J=4: 

3000.0 

4.31E-6 


TABLE(NI=2J4J=4): WAL2: 1=1,2 $alpha2 
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J=l: 

000.0 

1.00E-6 

J=2: 

1000.0 

2.26E-6 

J=3: 

2000.0 

3.28E-6 

J=4: 

3000.0 

4.31E-6 


TABLE(NI=2,NJ =4): WAL3: 1=1,2 $alpha3 


J=l: 

000.0 

1.00E-6 

J=2: 

1000.0 

2.26E-6 

J=3: 

2000.0 

3.28E-6 

J=4: 

3000.0 

4.31E-6 


T ABLE(NI=2,NJ=2): WYXX: 1=1,2 $ Yxx 

J=l: 0.0 1.0 

J=2: 30000.0 1.0 

TABLE(NI=2,NJ=2): WYYY: 1=1,2 $ Yyy 

J=l: 0.0 1.0 

J=2: 30000.0 1.0 

TABLE(NI=2,NJ=2): WYZZ: 1=1,2 $Yzz 

J=l: 0.0 1.0 

J=2: 30000.0 1.0 

TABLE(NI=2,NJ=2): WYXY: 1-1,2 $Yxy 

J=l: 0.0 1.0 

J=2: 30000.0 1.0 

TABLE(NI=2,NJ=2): WYXZ: 1=1,2 $ Yxz 

J=l: 0.0 1.0 

J=2: 30000.0 1.0 

TABLE(NI=2,NJ=2): WYYZ: 1=1,2 $ Yyz 

J=l: 0.0 1.0 

J=2: 30000.0 1.0 

$ 

$ 

$ Free up space 
!LINT = FREE ( ) 

ILTOT = FREE ( ) 

!L0 = FREE ( ) 

!L1 = FREE ( ) 

!L2 = FREE ( ) 

!L3 = FREE () 

!L4 = FREE ( ) 

!FLX = FREE ( ) 

!MTCC = FREE ( ) 

$ 

♦RETURN 

* ZDATi\ 

$ 

♦XQTAUS 

$ 

♦DC ALL (29 STRU ANAL) 

♦XQT EXIT 
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Table 1 - Temperature Dependent Properties Used in S81 and E43 Elements 
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Temperature E, psi v a, in/xn.°F 


Material 1 


0.0 

1.494 x 10? 

0.4000 

2.00 x 

10- 6 

1000.0 

1.300 x 107 

0.3000 

3.26 x 

10-6 

2000.0 

1.146 x 10 7 

0.2000 

5.28 x 

10-6 


Material 2 


0.0 

5.903 x 10 7 

0.1054 

1.00 x 

10-6 

1000.0 

3.712 x 10 7 

0.1034 

2.26 x 

10*6 

2000.0 

2.685 x 10 7 

0.1046 

3.28 x 

IQ 6 
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Table 2 - Comparison of Plate Stresses in the y Direction for S81 and E43 Elements 


Element 

a 

°S8!>P si 

b 

<*881’ P si 

<*E43» P si 

1 

-43.6 

-43.9 

-43.2 

2 

-40.2 

-40.6 

-40.2 

3 

-42.2 

-42.6 

-42.1 

4 

-45.1 

-45.4 

-45.0 

5 

-53.3 

- 53.7 

-53.0 

6 

-41.4 

-41.8 

-41.3 

7 

-42.0 

-42.4 

-41.8 

8 

-42.8 

-43.1 

-42.6 

9 

-47.1 

- 47.5 

-46.9 

10 

-51.1 

-51.5 

-51.0 

11 

-40.9 

-41.3 

-40.9 

12 

-42.1 

-42.4 

-41.7 

13 

-43.5 

-43.9 

-43.3 

14 

- 47.2 

-47.5 

-46.8 

15 

-50.7 

-51.1 

-50.7 

16 

-40.9 

-41.2 

-40.8 

17 

-41.0 

-41.4 

-40.8 

18 

-45.9 

-46.2 

- 45.7 

19 

-45.7 

-46.0 

-45.5 

20 

-50.9 

-51.4 

-50.8 

21 

-36.0 

-36.3 

-35.6 

22 

-43.9 

-44.3 

-43.9 

23 

-49.9 

-50.3 

-49.8 

24 

-48.5 

-48.9 

-48.5 

25 

-46.1 

-46.5 

-45.7 

26 

-29.1 

-29.3 

-28.6 

27 

-49.1 

-49.6 

-49.1 

28 

-53.1 

-53.6 

-53.2 

29 

-53.9 

-54.4 

-54.0 

30 

-39.2 

-39.4 

-38.6 


a Material properties interpolated over 1000 °F interval 
b Material properties interpolated over 100 °F interval 




Table 3 - Comparison of Analytical and EAL Displacements in a 1-D Rod 
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Position, in. 

EAL displacement, in. 

Eq. (4) displacement, in. 

% error 

0.00 

0.0 

0.0 

0.00 

0.05 

-0.231 x 10^ 

-0.241 x 10- 4 

4.15 

0.10 

-0.448 x 10- 4 

-0.464 x 10- 4 

3.45 

0.15 

-0.649 x 10- 4 

-0.668 x 10- 4 

2.84 

0.20 

-0.831 x 10- 4 

-0.851 x 10- 4 

2.35 

0.25 

-0.994 x 10- 4 

-0.101 x 10-3 

1.58 

0.30 

-0.114 x 10' 3 

-0.115 x 10-3 

0.87 

0.35 

-0.125 x 10- 3 

-0.127 x 10*3 

1.57 

0.40 

-0.135 x 10- 3 

-0.136 x 10*3 

0.74 

0.45 

-0.142 x 10-3 

-0.142 x 10-3 

0.00 

0.50 

-0.146 x 10-3 

-0.146 x 10-3 

0.00 

0.55 

-0.147 x 10-3 

-0.147 x 10-3 

0.00 

0.60 

-0.145 x 10-3 

-0.145 x 10-3 

0.00 

0.65 

-0.140 x 10-3 

-0.139 x 10-3 

0.72 

0.70 

-0.131 x 10-3 

-0.130x 10-3 

0.77 

0.75 

-0,119 x IQ' 3 

-0.118 x 10-3 

0.85 

0.80 

-0.104 x 10-3 

-0.102 x 10-3 

1.96 

0.85 

-0.839 x 10^ 

-0.826 x 10- 4 

1.57 

0.90 

-0.602 x 10- 4 

-0.591 x 10- 4 

1.86 

0.95 

-0.323 x 10- 4 

-0.317 x 10- 4 

1.89 

1.0 

0.0 

0.0 

0.00 


Table 4 - Comparison of Analytical and EAL Stresses in a 1-D Rod 


T 0 ,°F 

t l ,°f 

°EAL> 

°Eq. (6). P si 

% error 

200.0 

200.0 

-48.0 

-48.0 

0.0 

1000.0 

1000.0 

-2000.0 

-2000.0 

0.0 

200.0 

1000,0 

-619.8* 

b 

-647.3 

4.25 

200.0 

1000.0 

-646.7 

-647.3 

0.09 


a Material properties interpolated between endpoint temperatures 
b Material properties interpolated over 100 °F interval 
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Figure 1. - 



How chart of additional runstream needed to calculate temperature dependent 
properties. 
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and nodal temperatures. 


MSA. Report Documentation Page 

S r.KV AO">*V2* JIM' 


9. Performing Organization Name and Address 

Analytical Services and Materials, Inc. 

Hampton, VA 23666 

Old Dominion University Research Foundation ^ 

Norfol k , VA 23508 13. Type of Report and Period Covered 

12. Sponsoring Agency Name and Address 

National Aeronautics and Space Administration Contractor Report 

Langley Research Center in Sponsoring Agency code 

Hampton, VA 23665-5225 

15. Supplementary Notes 

D. E. Glass: Analytical Services and Materials, Inc., Hampton, Virginia. 

J. C. Robinson: Old Dominion University Research Foundation, Norfolk, Virginia 

C. J. Camarda: Langley Technical Monitor 


16. Abstract 

This paper presents a procedure to allow the use of temperature dependent 
mechanical properties in the Engineering Analysis Language (EAL) System for 
solid structural elements. This is accomplished by including a modular runstream 
in the main EAL runstream. The procedure is applicable for models with 
multiple materials and with anisotropic properties, and can easily be incorporated 
into an existing EAL runstream. The procedure (which is applicable for EAL 
elastic solid elements) is described in detail, followed by a description of 
the validation of the routine. A listing of the EAL runstream used to validate 
the procedure is included in the Appendix. 


18. Distribution Statement 

Unclassified - Unlimited 
Subject Category - 39 


page) 21. No. of pages 

39 


NASA FORM 1S26 OCT 86 


22. Price 

AO 3 


17. Key Words (Suggested by Author(s)) 

temperature dependent 
mechanical properties 
finite element analysis 


19. Security Classif. (of this report) 20. Security Classif. (of 

Unclassified Unclassified 


1. Report No. 2. Government Accession No. 

NASA CR-187470 

4. Title and Subtitle 

Procedure for Implementation of Temperature- Dependent 
Mechanical Property Capability in the Engineering 
Analysis Language (EAL) System 


7. Author(s) 

David E. Glass and James C. Robinson 











